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' Abstract We design a device that generates fields canceling out a known probing 

field inside a region to be cloaked while generating very small fields far away from 
the device. The fields we consider satisfy the Laplace equation, but the approach re- 
mains valid in the quasistatic regime in a homogeneous medium. We start by relating 
the problem of designing an exterior cloak in the quasistatic regime to the classic 
problem of approximating a harmonic function with harmonic polynomials. An ex- 
plicit polynomial solution to the problem was given earlier in [Phys. Rev. Lett. 103 

■ (2009), 073901]. Here we show convergence of the device field to the field needed 
to perfectly cloak an object. The convergence region limits the size of the cloaked 

■ region, and the size and position of the device. 

■ Keywords Cloaking • Laplace equation • Quasistatics • harmonic polynomial 
approximation. 
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1 Introduction 



> 
(N 

. Cloaking - preventing detection of objects from a probing field - has been the subject 

I of many recent studies, see e.g. the reviews IdtiZ]- A cloak can be active or passive 

. depending on whether active sources are needed to maintain the cloak. A cloak is 

said to be interior if it completely surrounds the object to be hidden and exterior 
otherwise. 

One approach to obtain passive interior cloaks is to exploit the invariance of the 
governing equations (e.g. Laplace, Helmholtz, Maxwell equations, . . .) to coordinate 
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transformations. This approach was introduced in ll6lEo [Jl4l[l5i [3i[7ll (see also refer- 
ences in |1, 7]) and is based on ideas first observed in Ahhough transformation 
based cloaking is set on solid mathematical grounds and has been demonstrated ex- 
perimentally in a variety of physical settings, the cloaks generated with this approach 
require materials with extreme properties that are usually approximated using spe- 
cially designed metamaterials. Unfortunately metamaterials used in electromagnetic 
transformation based cloaking are typically very dispersive, meaning that the cloak 
operates only in a narrow band of frequencies. Also losses in the cloak material gen- 
erate heat that can make the object detectable using infrared. Some recent results in 
generating broadband low-loss metamaterials have been obtained in |23]. In an effort 
to overcome the shortcomings of transformation based cloaks, various regularizations 
have been proposed (see [ 12] and references therein). 

Other passive interior cloaking methods include plasmonic cloaking (see 111] and 
references therein). Cloakin g m ethods that are passive and exterior include cloaking 
with complementary media II13I1 . cloaking by anomalous resonances 

EEH and 

plasmonic cloaking {72^. 

An example of an active interior cloak appears in |16] and uses sources contin- 
uously distributed over a closed surface surrounding the cloaked region in order to 
cancel out the incident field inside the cloaked region. 

Here we focus on an active exterior cloak for the 2D Laplace equation |8], which 
can be easily adapted to 2D quasistatics in a homogeneous medium. This scheme 
assumes the incident or probing field is known and uses one active source (cloaking 
device) to cancel the incident field in the cloaked region with no significant perturba- 
tion in the far field. Thus an object inside the cloaked region interacts very little with 
the probing field and becomes harder to detect. Active exterior cloaking has been ex- 
tended to the 2D Helmholtz equation in and to the 3D Helmholtz equation in 
1 1 in . Our approach assumes a homogeneous background medium and requires three 
(resp. four) devices or antennas to construct a cloak for the 2D (resp. 3D) Helmholtz 
equation. 

Our goal here is to rigorously justify the quasistatic cloaking method of |8]. Qua- 
sistatics refers to the Maxwell or Helmholtz equations in the long wavelength limit, 
where the governing equation is the Laplace equation. We start by describing the 
cloak setup in Section |2l Then in Section [3j we prove the existence of a solution for 
the 2D quasistatic active exterior cloak, based on a classic harmonic approximation 
result due to Walsh (see e.g. |5]). Unfortunately the existence proof is not construc- 
tive. We proposed a candidate constructive solution without proof in |8], supported 
by numerical experiments. In Section |4] we give the arguments behind this solution 
and prove that it does indeed solve the active exterior cloaking problem. 



2 Cloak setup and device requirements 

Three regions in are needed to describe our cloak setup: the region to be cloaked, 
the cloaking device, and the observation region. See Figure [T] (left) for an example 
setup. The main idea of our cloaking method is to cancel out an (assumed known) 
incident field mq inside the cloaked region while perturbing the far field only slightly. 
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Fig. 1 The effect of the inversion (Kelvin) transform w = 1 /z on the cloak geometry. The cloaked region is 
in red and the device sources ai'e all contained in the gray disk. The green region is the observation region, 
where the device field must be very small to avoid detection. 



Thus the total field inside the cloaked region is practically zero and the scattered field 
from any objects inside the cloaked region is reduced significantly. 

Here we consider the conductivity equation with conductivity one and a harmonic 
incident field uq (i.e. Auq = 0). Without loss of generality, we take as cloaked region 
the diskB(c,a) C M^, centered at c = (p,0) G M^, p>0, and with radius a > 0. As in 
js^, we consider one cloaking device located inside B{0,5), with 5 <C 1. The device 
generates a field u, harmonic outside B(0, 5). In order to cloak objects the device field 
u needs to satisfy the following requirements. 

1. The total field m + mq in the cloaked region B{c,a) is very small. 

2. The device field u is very small far away from the device, e.g. in the observation 
region \ B{0,R) for a large R>0. 

In order for the device to be exterior to the cloaked region, we must have 

p>a + 5. (1) 

Also the observation radius R needs to be large enough to contain both the device and 
the cloaked region: 

R>a + p. (2) 



3 Cloak existence 

The existence of a device field u having the desired cloaking properties to within a 
tolerance e is stated in the next theorem. 

Theorem 1 Let e > be an arbitrarily small parameter. Let also a > 0, c = (/?,0), 
p > Q and R satisfy the inequalities ([T]) and (|2]i. Then for a harmonic incident field 
Mq, there are functions g^:'^ and m : — > R such that 

Au^O, in r2\B(0,5), 
« = ^0, ondB{Q,5), (3) 
with \u\ < e in R^\B{0,R) and \u + uo\ < e in B{c,a). 
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Region 


z plane 


w = l/z plane 


Cloaking device 


B{0,S) 


]R2\B(0,1/5) 


Cloaked region 


B{c,a) 


B(c*,a) with c* = (/3,0), a = a/\p- - c?\ 






and/5 = /V(/-a^) 


Observation region 


R^\B{0,R) 


B(0, 1 /R) 



Table 1 The different regions in our cloak setup and how they are mapped by the inversion (Kelvin) 
transformation. 



The main idea of the proof of Theorem [T] is to relate active exterior cloaking to 
the problem of approximating harmonic functions with harmonic polynomials. We 
rely on the following classic result. 

Lemma 1 (Walsh, see e.g. 1 5], page 8) Let Kbea compact set in such that \ A' 
is connected. Then for each function w harmonic on an open set containing K and 
for any e > 0, there is a harmonic polynomial qfor which \w ^q\ < £ on K. 

We can now proceed with the proof of Theorem[T] 

Proof It is convenient to use complex numbers z—x + iy to represent points {x,y) £ 
S?. By applying the inversion (Kelvin) transformation w = l/z, the geometry of the 
problem transforms as in Table[T] (see also Figure[T]i. 

Thus the cloaking problem ^ is equivalent to finding functions go and u for 
which 

Au = 0, inB(0,l/5), 
u = go, ondB{Q,l/5), (4) 

with |m| < £ onB{0,l/R) and |m + mo| < £ onB(c*,a). 

Here £ is as in the statement of the theorem, go{z) = go{^/z) and the function uo{z) = 
Mo(l /z) is harmonic on \ {0}. 

Let Uq denote the analytic extension of mq in B(c* , a), obtained by adding / times 
its harmonic conjugate. Notice that since Uq is analytic, it can be arbitrarily well 
approximated by a polynomial, e.g. a truncation of the power series of Uq- Therefore, 
there is a polynomial Qq such that 

\Uo-Qo\<e/2, onB(c*,a). (5) 

For Uq this means that 

|mo-^o| < e/2, onB(c*,a), (6) 

where qo is the real part of Qq. Thus we may solve (|4| by first approximating the 
(inverted) incident field So by qo and then studying the following problem 

Au = 0, inB(0,l/5), 
u^go, ondB{0,l/5), (7) 

with |m| < £ onB(0, 1//?) and \u + qo\ <e/2 onB(c*,a). 
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After inversion, the conditions ([T]i and (|2]i necessary for having an exterior cloak 
become 

l/R< P - a, (the two disks B(0, 1 /R) and B{c* , a) do not touch), and 

(8) 

j3 + a < 1/5, (the two disks B(0, 1/5) and B{c*,a) do not touch). 
Therefore, there exists < <C 1 such that 



We can now apply Lemma[T|to the compact set K — B{0, l/R)UB{c* ,a) (which has 
a connected complement by virtue of ([8]l) and the function 



which is a harmonic function in the open set B(0, ^ + ^) UB(c*, a + (^) (a set contain- 
ing K). We obtain that there exists a harmonic polynomial q such that |q' — w| < e/2 
on K. A solution to ^ is then given by m = <7 and gQ=q on dB{0, 1 / 5). This implies 
the statement of the theorem. 

Remark 1 We assumed throughout this section that the incident field uq is harmonic 
on M}. This corresponds to a source located at infinity. Recall our method relies on 
approximating the Kelvin transformed analytic extension of the incident field Uq in- 
side the Kelvin transformed cloaked region B{c* ,a) by a polynomial Qo (see (fTTI) ). 
This approximation only requires analyticity of Uo inside the cloaked region B{c,a). 
Hence the results of this section and the construction of Section |4]below generalize 
easily to the case where the incident field uq is harmonic inside the observation region 
B{Q,R). This is the case where the sources generating the incident field are outside 
the observation region but not necessarily located at infinity. 

Remark 2 Clearly, Theorem[T]also holds when the device and cloaked region are not 
necessarily disks. The only requirements are that they be bounded, disjoint and that 
the complement of their union be connected (see Lemma[Tli. 

4 A constructive solution for active cloaking 

Although mathematically rigorous, the existence result of Theorem [T] does not give 
an explicit expression for the potential required at the active device (antenna). To give 
an explicit harmonic solution to problem (O, we first simplify the problem in Sec- 
tion |4T| Then we give a candidate solution to the simplified problem in Section l4!2l 
in the form of a Lagrange interpolation polynomial. A better solution is constructed in 
Section l43] bv averaging several Lagrange interpolation polynomials. The resulting 
polynomial turns out to be a Hermite interpolation polynomial. Then in Section l4~4l 
we show that this Hermite interpolation polynomial solves (|4]i (and thus the cloaking 
problem (O) provided its degree is sufficiently large. This convergence study reveals 
constraints on the size of the cloaked region and the device that are due to the partic- 
ular solution we construct. 



(9) 
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4. 1 Simplifying the problem 

In the proof of Theorem [1] we related the cloaking problem Q to the problem of 
approximating a polynomial Qq with an analytic function V such that for some e > 0, 

|y| < e in 5(0,1/7?) and |y + eo| <einB(c*,a). (11) 
Now consider the problem of finding an analytic function W such that for some e' > 0, 

|1 -W| < e' mB{Q,l/R) and |W| < e' inB(c*,a). (12) 
Assuming we can find an approximant W in ( fT2] i with e' = e/M and 

M= sup iGoWl, (13) 

zGB(c*,a)UB(0,l/ff) 

a solution to (fTTl i is then V = — 2o(l— W'), which is analytic because the product of 
two analytic functions is analytic. 

For illustration purposes we fast forward to Figure[3] where we give an example of 
a function W with the approximation properties ( fT2] i. The function W is a polynomial 
whose motivation, derivation and analysis are the subject of the remainder of this 
section. 

In order to use such a function W for cloaking, assume Qq{\/z) is the harmonic 
incident field. Then the device field needed for solving the cloaking problem (O is 
the real part of the function U{\/z) = —Qo{^/z){l — W{1 /z)) (after having undone 
the Kelvin transformation we used for the analysis). The actual device field is illus- 
trated in Figure|4] On the left, a scatterer perturbs the incident field and can be easily 
detected. On the right, the device field (based on the function W of Figure [3] is acti- 
vated and suppresses the incident field inside the cloaked region, making the object 
undetectable for all practical purposes. 

4.2 A first candidate polynomial from Lagrange interpolation 




Fig. 2 Left: sample interpolation points for the interpolation polynomial with ;i = 5, = and 
l/A = n/'i. Right: the modulus of the polynomial Pi/i.ij/ with n = \0, (^i = — = n/W and /3 = 4. The color 
scale is logarithmic and the interpolation nodes are indicated by the interpolation values. 

We present a polynomial solution to (fT2l i based on Lagrange interpolation. This is 
an intermediary step to motivate the explicit solution to (fTST i given later in Section l43] 
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The idea applies only to the case where a=R = I and /3 = p/{p^ — a^) > 2. The can- 
didate solution is a polynomial that is one at « equally distributed points on dB{Q, 1 ) 
and zero at n equally distributed points on dB{c* , 1). The motivation being that by 
surrounding both and c* = (j3 , 0) by n points where the polynomial has the desired 
values, we hope to get close to a polynomial satisfying (fT2] l. 

To be more precise, let us introduce the following family of 2n nodes {e"^Wj , j3 + 
e'^Wj}"!^- Here and Xj/ are two arbitrary angles and wj ~ exp[2i7tj/n], for j = 
0, . . . ,n — 1. Define the polynomial p^^ as the unique polynomial of degree 2n — 1 
satisfying. 



p<^^^,{e'^wj) = 1 and p^.^(j3 +e'^Wj) = 0, for j = 0, . . . ,n - 1. 



(14) 



An example of the interpolation nodes and the values of p,j,,if, is shown in Fig- 
ure lltleft). 

The polynomial p,;,.^ is unique and can be written explicitly as 



n-1 
m=0 



(15) 



where q^^y/^miz) are Lagrange interpolation polynomials (see e.g. 02411 ) defined for 
m = 0, ...,«— 1 by 



n-1 

n 



jJ)e'^w,„-(j3+e'Vw;) 



(16) 



1=0, j^m J 

or alternatively by their interpolation properties 

q(l>,\tr,m{e'^Wj) ^ 8mj, and ^0^^^„,(j3 + e'"''wy) = 0, for ; = 0, . . . ,n - 1. (17) 

Here 5mj = 1 if m = j and otherwise is the Kronecker delta. Straightforward calcu- 
lations give the expression 



,\l/,m{z) 



iz-py-e 



Z — < 



1 



n{w,ne"l>) 



i6 \n— 1 



(18) 



which will be used later in Section l431 

We state the following symmetry property of the polynomial for later use. 

Lemma 2 For any angles <j) and the polynomial p^,\ji has the following symmetry 
property: 

P^,v{z)+p^,+Ti,^+n{^-z) = l. (19) 



Proof Equation ([19]) follows from noticing that for y = 0, . . . , n — 1, 

PV+K4+AP - iP+e'^^^j)) = Pv+^4+^i^'^'^^'''^^j) = 0, and 

P\li+7c,(li+7iiP ~ e"^Wj) — Pii,+n.ip+7ciP - 



J{0 + 7i] 



(20) 



Hence the polynomial p<j),\ir{z) + P\ji+K,^+K{P — z) — 1 must be identically zero be- 



n-l 



cause it is of degree 2« — 1 and has 2n roots {e"^Wj,^ + e'^Wj}"^ 
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An actual polynomial p^^f, is shown in Figure fright). Unfortunately this poly- 
nomial is not a good solution for problem (fTSl i as the regions where p,j, ^, w 1 and 
Ptl),\l/ ~ to within a certain tolerance (say 1%) are relatively small. Changing (j) and 
xj/ does not give a significant improvement. However these polynomials are the build- 
ing block for the ensemble average polynomial solving (fT2] i that we present next. 



4.3 The ensemble average polynomial 





Fig. 3 The modulus of the ensemble average polynomial (p) {z) for n = 12 and j} = I. The device field 
used for cloaking is (p) (1/z). Within 1% accuracy, the polynomial (p) is close to one inside the dashed 
white circle and close to zero inside the solid white circle. The boundary of the convergence region Dp of 
(p) as n — > oo is the peanut shaped curve in red (see Theorem|2j- The color scale is logarithmic from 0.01 
(dark blue) to 100 (dark red), with light green representing 1. 



In an effort to obtain a polynomial solution to problem (fT2l i we calculate the 
ensemble average of the polynomials with respect to the two phase factors 
0,1^ e [0,271], that is 



In r2n 



''^^ " (2^ Jo Jo ^^■v(^)«''?'^V^- (21) 

We prove in Theorem|2l using the next lemma, that indeed (p) is a solution for (fT2l) . 
An example of such polynomial for j3 = 1 and « = 12 is given in Figure[3] 

Lemma 3 The ensemble average polynomial defined in (I211 i has the expression 

<.>.)^(.-f)'i:a)'^. 

Proof We first use the Cauchy residue theorem to compute the integral 

1 f^^ jz-pf-e'V" __L f {z-pT-w" dw 



2;rJo {e'^w„,-^Y~e't"' ^ lin J\n-\=i {e"t'w„r - P)" - w" w 
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Fig. 4 Real part of the total field with the cloaking device active (right) and inactive (left), for an incident 
field «o(z) = z and n = 12. The solid white, dashed white and red curves are the Kelvin transforms of 
the respective curves in Figure[3] The solid black disk is an almost resonant scatterer with radius r = 0.1, 
located at z = 1.05 and with dielectric constant e = — 1 + 10^^, chosen to be plasmonic with a negative 
value close to —1 to amplify its effect. The solid black curve is the contour \u\ = 100. The color scale is 
linear from -10 (dark blue) to 10 (dark red). 



since the integrand has a single simple pole at w — in the disk |w| < 1. Then by 
plugging (l23T l into the expression for q^,\j/,m we get that 



1 f^^ 



(z-j3)" z"-< 



1 



2nJo -'y^-^'-- (e'^w,„-j3)«z-e'^w„,n(w„e'^)«-l' 
Recalling that pij, ^, is the sum (fTSl ) of qii).\f/.m we can write 



{P)iz) 



1 



Now all n terms in the previous sum are identical, therefore 

1 f^"" {z-PT zf-e'^" 1 



{P)iz) 



2k Jq {e'^-^)" z-e'<l> e'<l'(«-i) 
_L f (- BY 

2iKj\„\=i ' (z-w)(w-j3)« w" 



2iK 

(z-nr 

2iK 



^-1 
=1 \ w" 

n-l 



1 



(z-w)(w-j3)« 



dw 



w|=ij^)W^-+i(w-j3)« 
1 



dw 



(w-j3)« 



n.-=0 



(24) 



(25) 



(26) 



where we used Cauchy's theorem in the last equality of ( |26] |. The desired expression 
(l22T i follows by straightforward algebraic manipulations of ( 
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Remarks By using elementary algebraic manipulations and ( |26] |. it is possible to 
show that (p) is the Hermite interpolation polynomial 1124] of degree 2« — 1 that is 
uniquely defined by the 2n interpolation conditions 

(;,)(0) = 1, W(/3)=0, and (;,)(^) (0) = (jS) = 0, for j = 1, ...,«- 1. (27) 

Notice that the ensemble average polynomial inherits the symmetry property ( fT9t 
for that is 

{p){z) + {p){p-z}^l. (28) 

This symmetry property means that by design, the polynomial gives as good an ap- 
proximation to one near the origin as the approximation to zero near j3 . 



4.4 Asymptotics of the ensemble average polynomial 

We now study the behavior of the polynomial (p) (defined in (|2TI) ) as n °o. The 
following result shows that the polynomial (p) solves the problem (fT2l i. and gives 
limits to the size of the cloaked region. 

Theorem 2 The ensemble average polynomial (p) can be written as 

<">4^i:ii(f(-t))'a-f)- - 

The polynomial (p) (z) converges asn if and only if z belongs to the convergence 
region 



D/3 = |zeC, |z^-j3z|<^^|. (30) 
The convergence is uniform on compact subsets of Dp to the function 

= "f'^^^"- (3.) 
I U otherwise. 

For large enough n, the polynomial (p) solves (fT2l i if and only if 
1 B B 

-< — ^ and a< — ^ . (32) 

R 2V2 + 2 2V2 + 2 

Proof Consider the function 

fn{t) = {i-tr"ft^(''+-'r^) (33) 
7=0 \ J / 

where for any positive integers m and p, (™) = jTi0rjPj[ - Note that, from (l22l) we have 
/„(f) = (/?)(/3f), forallf eC. (34) 
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Then for all f 7^ 1 we obtain, 

fn+i{t) 



(1-0 



n+1 



)=0 



n + ] 



.7=1 



n + j-l 



n + j-\ 



" ,./«+/■- 1\ r " ifn + j-1 



n-l 



k=Q 



n + k 
k 



fn+ljt) 



l-t 



2« 



2«-l 
« 



In the above equation we used the recurrence relation 



m — 1 



for any integers m,p > 0. 



From ( l35T l, for any integer n > 1 we obtain, 

/„+i(0-/„(f)-(i-0"f"+' 

= /„(r)-(l-0"f" 



+ (l-f)"f" 
1 



2n-l 
n 



/ - - ) , for all f 7^ 1 . 



In ( |36] | we used the identity 



2n- 1 
n 



, for every integer n > 1 . 



From the first order Unear recurrence (l36b we obtain. 



n-l 



/t=0 



(35) 



(36) 



(37) 



and this is valid for all n > 1 and all f e C (as ( |37] ) which was initially obtained 
for t I checks also for t — 1). The final expression ( |29] | follows from substituting 
t = z/15 in ^ and using ( l34l l. 

Notice that the polynomial (p) is in fact the n-th order partial sum of the following 
infinite sum. 



k=0 



2k' 



1 



By the ratio test this series converges uniformly on compacts subsets of the region 
Dp (defined at ( l30b ) to a limit function (p and diverges in C\D^. From the uniform 
convergence of (p) we deduce the analyticity of (p in and by using the Taylor 
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expansion around the origin for (p, the Remark[3] and the symmetry property (l28T l we 
obtain convergence to the function ( |3TI ) inside Dj^ . 

We now study the convergence region Dp in order to show that the constraints 
( l32l l are necessary and sufficient for (p) to solve (fT2l l. First notice that the definition 
of the region Dp and simple algebra reveal that 



(38) 



4 < ^ 4^"" e (£)« n {z e C,2;;e(z) < /3}), and 

2V2 + 2 R ^ P ^ 

a < — ^ j3 + a e (Dp n {z e C,2/?e(z) > /3}). (39) 
2v 2 + 2 



Next we show that 

^ < —I ^ B(0, 1//?) (E {Dp n {z e C, 2/?e(z) < jS}), and (40) 

2v2 + 2 

a < —I — ^ B(c*,a) (E {DpCiizeC, 2Re{z) > p}), (41) 
2v 2 + 2 

where (g is the classical symbol for compact inclusions. By using the equivalences 
(l38T l and i39[ it is easy to check that the inclusions in ( l40l i and (HTt imply the con- 
straints (|32] ). To show the implication (=>), we first show that for any two positive 
real numbers l,q with 2max{l,q} < /3, we have 



le-"' eDp-^ B{0,1) <£ {Dpn{z e C,2Re{z) < p}) and (42) 
P+qeDp ^B{c*,q) (E (D|3n{zeC,27;e(z) >J3}). (43) 

Let us first show the equivalence (l42l l. The sufficiency (<^=) is immediate. For the other 
implication (=>), we can use the definition of Dp to show that for any 9 E [— tt, n] we 
have. 



fi2 



4 

^ (/V-j30(/V'''-j30<T7 

16 



R4 

<^ Z'^ + i3¥-2Z^/3cos0-^ <0. (44) 

16 

Since we assumed le^'" G D^, equation ( l44l i immediately impUes that 

l'^ + l5^f + 2Pl5-i—<0. (45) 
16 

Consider the even function / : [— ;r, ;r] — > M defined by, 

fi4 

= /4 + j3¥-2/^j3cos0-^. (46) 

16 

Observe now that, because I > 0, its derivative f'{9) — 2/^/3 sin 9 has the signs, 

>0 fore e [0,;^] and < for e [-;^,0]. (47) 
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Note that from inequality ( l45l l and the definition ( |46] | of /(0) one immediately ob- 
tains 

f{~7t)^fi7Z)<0. (48) 

Then the signs of f'{9) in (|47] | together with the particular values of f{9) in ( |48] ) 
imply 

< max{/(7r),/(-7r)} < 0, for all 9 e [-7r,7r]. (49) 
Because of equivalence (|44] | we conclude from inequality (|49] l that 

Ze'" e D|3, for any e [-:t,7i;]. (50) 

From the conditions on /, we have that 2Re[le'^) < 2/ < j3 and by using this in dSOl ) 
we obtain 

le''^ e {Dp n {z e C,2;;e(z) < /3}), for any 9 e[-n,n]. (51) 

Inclusion ( fSTT l together with the convexity of n {z G C,2/?e(z) < /3} implies that 

5/(0) (E D|3 n {z e C,27;e(z) < j3}. 

This establishes the equivalence (l42l l. From the definition of the set D^, by simple 
algebraic manipulation we obtain 

P+qGDp^ qe-''' € Dp . (52) 

Equivalence ( |52] | clearly implies that ( |43] | follows from (l42l) applied to instead of Z. 
Finally, observing the fact that the constraints ( |32] | imply 2max{ ^, a} < /3 and using 
equivalences dSST l. ( [39] l. ( |42] | and (l43T l for ^ and a instead of / and q respectively, we 
obtain the desired equivalences ( |40] | and (l4Tl i. By using the uniform convergence of 
the polynomial (p) to the function x {z) in Z)j3 , and equivalences ( l40b and (|4TI) we 
obtain that the constraints (l32b are indeed necessary and sufficient for convergence of 
ip)- 

Remark 4 The expression ( |29] | of the ensemble average polynomial could also be 
obtained by generalizing to distributions a theorem by Ramharter f^T] (which is in 
turn a generalization of a result due to Berger and Tasche [2J). To remain concise, we 
prefer to include a direct proof. 

5 Summary 

For the Laplace equation we have shown the existence of a device capable of cloaking 
a region exterior to the device, assuming a priori knowledge of the incident field. 
The proof relies on a non-constructive harmonic function approximation result. The 
theory does not constrain the size and relative positions of the device and cloaked 
region, as long as they are bounded, disjoint and the complement of their union is 
connected. Although the construction of such a cloaking device is clearly not unique, 
we presented earlier in |Q] a construction based on an explicit polynomial. Here we 
rigorously justify this construction and show that the constraints ( [32b must be satisfied 
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in order to have a proper active exterior cloak. Because of the constraints ( |32] ). the 
current strategy fails to cloak large objects (a large) unless they are sufficiently far 
from the origin (j3 large enough). In |11] (see Conjecture 1), we present without 
proof, as a conjecture, an extension of Theorem |2] which gives a wider choice of 
cloaks and that is supported by numerical experiments. 
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